Time-evolution of tripartite quantum discord and entanglement under local and 

non-local random telegraph noise 
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Few studies explored the dynamics of non-classical correlations besides entanglement in open 
multipartite quantum systems. Here, we address the time-evolution of quantum discord and en- 
tanglement in a model of three non-interacting qubits subject to a classical random telegraph noise 
in common and separated environments. Two initial entangled states of the system are examined, 
namely the GHZ- and W-type states. The dynamics of quantum correlations results to be strongly 
affected by the input configuration of the qubits, the type of the system-environment interaction, 
and the memory properties of the environmental noise. When the qubits are non-locally coupled to 
the random telegraph noise, the GHZ-type states partially preserve, at long times, both discord and 
entanglement, regardless the correlation time of the environmental noise. The survived entangled 
states turn out to be also detectable by means of suitable entanglement witnesses. On the other 
hand, in the same conditions, the decohering effects suppress all the quantum correlation of the W- 
type states which are thus less robust than the GHZ-type ones. The long-time survival of tripartite 
discord and entanglement opens interesting perspectives in the use of multipartite entangled states 
for practical applications in quantum information science. 



I. INTRODUCTION 

Entanglement is a peculiar feature of quantum me- 
chanics which constitutes a valuable resource for a num- 
ber of features of quantum information processing, such 
as quantum communication, cryptography, and compu- 
tation However, the major limitation of the use of 
entanglement in practical applications is due to the un- 
avoidable interaction of the real quantum systems with 
their surroundings, resulting in decoherence processes, 
and, as a consequence, in degradation of the entangled 
states 0. Therefore, many research efforts focused on 
the investigation of the decoherence processes and of the 
dynamics of entanglement in different quantum systems 
ranging from quantum optics to nanophysics Q. 

The decohering effects of the environments with differ- 
ent memory properties have widely been studied in bi- 
partite systems where entanglement can exhibit peculiar 
features including sudden death and transitions, revivals, 
and trapping 04a, • On the other hand, only recently 
the time evolution of the entanglement has begun to be 
analyzed in multipartite open quantum systems, since 
in such systems the characterization and quantification 
of the quantum correlations still represent a hard task. 
Specifically, different approaches have been used [f| ®- 
Il3| . Some of them involve the generalization of entan- 
glement measures ranging from negativity @ to concur- 
rence [f| |T(| or Bell nonlocality disequalities 11 1 com- 
monly adopted in bipartite systems. Multipartite entan- 
glement has also been investigated by means of suitable 
entanglement witnesses (EWs), namely observables that 
can detect the presence of the entanglement itself [T243] ■ 
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Apart from entanglement, quantum states can exhibit 
other correlations not present in classical systems. In 
particular, in the last years quantum discord (QD) has 
received an increasing attention as a mean to estimate 
the quantum correlations encoded in a system and po- 
tentially usable for applications in quantum informa- 
tions processing }15H19| . Indeed, it was shown theoret- 
ically [2(j, and later experimentally [2lj], that bipartite 
separable mixed states with nonzero discord may pro- 
vide computational speedup compared to classical states 
in some quantum computation protocols. Under suit- 
able conditions quantum discord is also found to be more 
robust than entanglement in noisy environments [l7| . 
While the characterization of quantum discord in terms 
of the discrepancy between two quantum analogues of the 
classical mutual information has well been grasped and 
widely used in bipartite systems, its extensions in multi- 
partite systems is still discussed and tackled with differ- 
ent approaches [22j-|27j. Recently, Giorgi et al. [26| intro- 
duced a measure of genuine total, classical, and quan- 
tum correlations relying on the use of relative entropy 
to quantify the distance between two density matrices. 
Like oth er p roposed measures of correlations in quantum 
systems [Ij, 0, H3| , the approach of Ref . (2(| involves dif- 
ficult extremization procedures over operators or states 
which makes the calculations very hard. This justifies 
the scarse number of works exploring the time evolution 
of quantum correlations in multipartite open systems. 

In this paper, we investigate the dynamics of the quan- 
tum correlations in a physical model consisting of three 
qubits, not interacting among each other, and subject to 
a classical random telegraph noise (RTN). Two different 
initial configurations are examined, namely GHZ and W 
Werner type-states. The influence of the classical envi- 
ronmental noise on the system is described by means of 
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a stochastic Hamiltonian with a coupling term mimick- 
ing a random telegraph signal. The dynamics of the two 
qubits is evaluated by averaging the time-evolved states 
over the noise. Even if in similar systems of two qubits 
the time evolution of quantum correlations has already 
been investigated [28- 31], the analysis of the tripartite 
regime seems quite promising since here the correlations 
derive from a more complex scenario, and, as a conse- 
quence, more peculiar phenomena are expected. 

The use of a stochastic Hamiltonian allows us to eval- 
uate the dynamics of qubits in both separated common 
environments in a straightforward way. While in the for- 
mer case each subsystem is locally coupled to a source of 
RTN, in the latter a non-local interaction between the en- 
tire system and the environment occurs, being the qubits 
collectively coupled to the same source of RTN. This gives 
us the chance to analyze how the local or non-local na- 
ture of the qubits-environment coupling and the memory 
properties of the environment itself affect the dynamics 
of quantum correlation and its possible sudden death, re- 
vival, or trapping. From this point of view, the choice of 
a classical environment does not represent a strict con- 
straint of the model, since a quantum environment not 
affected by the system or influenced in a way that does 
not result in back-action, can be mimicked by a classical 
noise [32[ that was already shown to be able to lead, in 
bipartite systems, to sudden death and revival of entan- 
glement and discord [H, |3l|, HH . 

Here, we focus on the quantum correlations in the form 
of tripartite entanglement and quantum discord. Being 
the characterization of tripartite mixed entangled states 
still debated, two different estimators are adopted to eval- 
uate analytically the entanglement in agreement with 
what was done in other works [H, [H| . One estimator 
is the tripartite negativity (3~3 | , which was shown to play 
a key role in quantum information protocols since states 
exhibiting nonzero tripartite negativity are distillable to 
GHZ states. The other estimator relies on the concept 
of detectable entanglement and uses suitable EWs (14j . 
The quantum discord is numerically quantified by using 
the approach developed in Ref . [26[ . Unlike the tripartite 
negativity and the expectation value of the EWs, in this 
case we need to use numerical techniques apt to optimize, 
over positive operator-valued measures (POVM's) and 
subsystems, the conditional entropies appearing in the 
definition of discord. The comparison between the time 
evolution of these estimators for two different initial con- 
figurations of the qubits can represent a valid guideline 
to shed some light on the discrepancies and the analogies 
between the two forms of correlation in our model. 

The paper is organized as follows: in Sec. |Hl we intro- 
duce the estimators used to quantify tripartite quantum 
discord, entanglement, and detectable entanglement. In 
Sec. IIII1 we illustrate the physical model consisting of 
three qubits subject to a classical RTN within common 
or different environments. Sec. lIVI reports the time evolu- 
tion of quantum correlations for two initial configurations 
of the qubits in the case of local and non-local system- 



environment interaction. Finally, conclusions and discus- 
sions are summarized in Sec. |Vj 



II. QUANTUM CORRELATION MEASURES IN 
TRIPARTITE SYSTEMS 

In this section, we briefly expose the correlations mea- 
sures adopted in our work to quantify and characterize 
the tripartite quantum discord and entanglement. 



A. Tripartite Quantum Discord 

Here, we illustrate the criterion adopted to quantify 
the genuine total, classical, and quantum correlations 
in a tripartite system. By following the approach used 
in Ref. [2g|, the genuine total correlations T^{p), of a 
mixed state p can be expressed as: 



T^(p)^T( p )-T {2 Hp) 1 



(1) 



where T(p) is the quantum extension of the Shannon 
classical mutual information 



T(p) = S(p A ) + S( PB ) + S(p c ) - S(p) 



(2) 



with S(pk)=— Trfp/f hi2 px) indicating the the von Neu- 
mann entropy of px ■ T' 2 -* (p) is given by the maximum of 
the total bipartite correlations among two qubits of the 
system, namely 

T (2) (p) = max [lW(p AtB ),l i2) (PA,c)a i2 HpB,c)} (3) 

with l^\px. J )=S{p I ) + S(pj) - S(pu) denoting the bi- 
partite quantum mutual information of the subsystem 
consisting of the two qubits I and J. By inserting 
Eqs. © and (J3J) in JTJ) , one finds 

t^(p) = 

l^(pA,Bc)^ {2) (pB,Ac)a i2) (PCAB)} , (4) 

that is the total correlations of the state p are the low- 
est bipartite mutual information between the one-qubit 
part and the two-qubit part. As indicated in Ref. [26j j . 
T^\p) estimates the amount of genuine tripartite cor- 
relations that cannot be accounted for considering any of 
the possible subsystems. 

Analogously with T^(p) , also the genuine classical 
correlations J 1 ^ (p) can be defined in terms of the differ- 
ence between the total classical ones J{p) and the max- 
imum among bipartite correlations J^ 2 '{p): 



J {3) (p) = J(p)-J {2) (p)- 



Here 



J{p) 



max \S(pr) 
i,j,k L 



S(pi\j) + S{p K ) - S{p K \ij) 



(5) 



(6) 
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denotes the maximum over the six possible quantum ex- imized on the complete measurement on two subparties 
tensions of the conditional entropy which, in turn, is max- where 

I 



Prn 



s (Pi\j) = min y2PmS(pi\El) with p m = Tr/j (E^pu) ,Pi\ E J 
{Ej} m 

S(pK\u)= min y^Pi m S{p K \ E .i E i) with p lm = Tr IJK (EjE^p) , P k \e£e' 



TruEfEjp 

Pirn 



(7) 



In the above expressions, {Ej}, and {E^} indicate the 
POVM's acting on the qubits / and J, respectively. On 
the other hand, J^^(p) corresponds to the maximum of 
the classical bipartite correlations among two qubits of 
the system 

J {2) (P) = 

max \jW (p A}B ),jW ( PAC ),JP) (p fl c )j , ( 8 ) 



with 



where ( W ,j)=max [j£}{pi,j)> Jj 2 }(pi,j) 

J?]{PI,J)=S{ PI )~S{ PIV1 ). 

In agreement with the usual definition for two-qubit 
systems [HI, HU, the tripartite quantum discord V^\p) 
can be expressed as 



v^( p ) = t (3 Hp)-J [3) (p)- 



(9) 



Given the optimization procedures over the indices and 
POVMs appearing in the formal expression of T^ z \p) 
and J^{p), the evaluation of Eq. © is not immedi- 
ate. For symmetrical states, namely those systems whose 
state is invariant under the permutations of the three par- 
ties, it can be shown that the tripartite total, classical, 
and quantum correlations can be evaluated by means of 
simple expressions. After a straightforward calculation, 
one finds 



r^(p)^s(pc) + s(p AB )-s(p) 

J {3) (p) = S(pc)-S(p clAB ) 

V^{p) = S( PAB ) + S( Pc \ab) - S(p). 



(10) 



These are formally equivalent to the forms reported in 
Ref. [27| but here, in agreement with the approach devel- 
oped by Giorgi et al. [26| the von Neumann entropy of the 
measurement-based conditional density operator pc\AB 
is minimized over the set of two-qubit POVM's factoriz- 
able in terms of single-qubit POVM's E^(g>E^. The use 
of the latter as optimal POVM in the definition of clas- 
sical tripartite correlations has been discussed [27|- In 
our numerical approach the von Neumann entropy of the 
measurement-based conditional density operator given in 
Eqs. ([7]) has been estimated by using projective measure- 
ments described by the complete set of orthogonal projec- 



I 

to, n = 0, 1) defined by the orthogonal states 

|CU ( b)) = cos04 (s) |O) + e l *^> sin0 A(B) |l) 

|1a<b)> = e"^< B » sin0 A(B) |O) - cos^ (B) |l>. (11) 

The conditional entropy S(pc\ab) has been minimized 
over the phases 6a, Ob, ®a, and & B . 

In the followings, we will use Eqs. ([TU)) to evaluate the 
dynamics of quantum discord since at any time t the 
system investigated is invariant under the permutations 
of the three qubits. 



B. Tripartite Negativity 

In order to quantify the three-qubit entanglement of 
a mixed state p, we use the tripartite negativity Af^ 
which is given by |34j |: 



AA (3) (p) - VMa-bcMb-acM 



C-AB- 



(12) 



tors II 



ABi 



,0B,$A,$B)=Y^mn\ m An-B}(nB'mA\ (with 



It is the geometric mean of the bipartite negativi- 
ties Mi-jk (where I=A,B,C and JK=BC, AC, AB), 
which, in turn, are defined as Afi-jK=Y^i \ a i{p T ')\ ~ 1j 
with cii (p Tl ) the eigenvalues of the partial transpose p Tl 
of the total density matrix with respect to the subsystem 
/. For symmetrical tripartite quantum systems, Af^ re- 
duces to the bipartite negativity of any bipartition of the 
system. Recently, Af^ 3 ' has also been used to estimate 
the entanglement of tripartite bosonic and fcrmionic sys- 
tems (35| . 

Even if the positivity of the tripartite negativity en- 
sures that the state of the system under investigation is 
not separable and distillable to a GHZ state, Af^ cannot 
be used to classify and fully characterize the entangle- 
ment of general mixed tripartite states 3J] Apart from 
pure states, null tripartite negativity could indeed not 
imply the absence of entanglement. 



C. Entanglement Witnesses 

Multipartite entanglement can also be investigated by 
means of entanglement witnesses (EWs), namely suit- 
able observables that, at least in principle, can be exper- 
imentally implemented to detect the presence of entan- 
glement [l2T - fl4j . Following the classification of Ref. [lij , 
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four classes of three-qubit mixed states can be defined: 
separable states (S), biseparable states (B), W states, and 
GHZ states. Note that each class includes within it the 
previous classes as special cases. EWs are suitable ob- 
servables whose expectation value is zero or positive for 
all the states belonging to a specific class and negative 
for at least one state of higher class, namely a more inclu- 
sive class. Thus the use of EWs allows one to determine 
which class a given state belongs. 

Specifically, in this work we will use those EWs that 
permit to identify whether a state is in the W-B class, 
namely a state with true tripartite entanglement either of 
the GHZ-type or W-type and not biseparable. As shown 
in Refs. [12, [HI, GZH-type states in the W-B class can 
be detected by means of 



Hamiltonian 131] 



%2 = -I-Pghz, 



(13) 



where Pghz is the projector onto \GHZ). On the other 
hand, for W-type states, we will use the EW 



in 



-I - A, 



(14) 



with Pw the projector onto \W). Negative expectation 
values of the observables of the Eqs. (flU)) and (fl"4"|) in- 
dicate the appearance of tripartite entanglement exper- 
imentally detectable in the system. However, zero or 
positive values do not guarantee the absence of entangle- 
ment. 



III. THE MODEL 

Here, we describe a model consisting of three non- 
interacting qubits subject to an environmental classical 
RTN. The system-environment coupling is here analyzed 
in two different conditions. In the first one, each qubit 
locally interacts with its environment. In the other case, 
all the qubits are coupled with a unique common source 
of RTN thus mimicking a non-local interaction between 
the qubits themselves and environment. In both con- 
figurations, the dynamics of the system is ruled by the 
Hamiltonian 

H(t) = H A {t)® Ibc + n B {t)® I ac + H c {t)® Iab, (15) 

where Ijk is the identity operator in the subspace of the 
two qubits J and K and %z,(t) denotes the single-qubit 

I 



H L (t) = eI L + vq L {t)a xl 



(16) 



with a X L indicating the Pauli matrix of the subspace of 
the qubit L. e is the qubit energy in the absence of 
noise (energy degeneracy is assumed), v is the system- 
environment coupling constant. The RTN is here intro- 
duced by means the stochastic process denoted by Tji,(t) 
describing a fluctuator randomly flipping between the 
values -1 and 1 at rate 7. The two-particle form of the 
model defined by Eqs. (fT5|) and (|T6| has already been 
used to evaluate the time behavior of entanglement and 
quantum discord [3l| . Furthermore its extension to qu- 
dits case allows also to analyze the quantum walks of 
not-interacting particles in one-dimensional lattices [36| . 

Given the random nature of the telegraph processes, 
the time-dependent Hamiltonian of Eq. (|15l) leads to a 
stochastic evolution of the quantum states. In order to 
obtain the time-evolution of the system under the in- 
fluence of RTN, we need to average over different noise 
configurations [30l. l3l|. Therefore, the dynamics of the 
system density matrix can be written as 

p(t) = (U({ V },t)p(0)UH{ V },t)) {ri} (17) 

where p(0) is the initial state of the three qubits and 
U ({r]},t)=UA(riA,t) <E> U B (r]B,t) ® U c {f]c,t) indicates 
the unitary time evolution operator of the system at time 
t for a given noise configuration {77}— {77^4, 773,770} (where 
T]A{t) 7?s(t) 7^ T)c(t) in the case of qubits locally cou- 
pled to their environments and r)A(t)=r)B{t)=T)c(t) for 
non-local interaction of three subsystems with a com- 
mon source of noise). The single-qubit time-evolution 
operator Ul (VL,t) can be written as: 



Ul (VL,t) = exp 



-1 1 n L (t')dt' 




exp (— iet) 



cos 4>L{t) i sin (j>L(t) 
ism(pL(t) cos 4>L(t) 



(18) 



where h=l and 4>L{t)—-v f Q T]L(t')dt' is the random phase 
picked up during the time interval [0,i\. The evalua- 
tion of the time-evolved density matrix of Eq. (|17l) re- 
quires the estimate of the averaged terms of the type 
([cos(/>(/j)] m [sin0( y ;)]' c ). These can be expressed in terms 



of the average of the phase factor 
explicit form reads [37j : 



(e m< ^*'), whose the 



(t) ) = (cosn^(t)) 



+ ilsinn 



(19) 



with 



(cos n#(t)) = G n (t) 
(sin n0(f)) = 0, 



e 7t cosh (S ni ,t) + j 2 - sinh {8 nv t) 



cos (5 nv t) + sin (8 nv t) 



for 7 > nv 
for 7 < nv 



(20) 
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where 5 nu ^yJ\"y 2 — (m>) 2 \ (with m,k,n € N). Here, we 
examine the time-evolution of two classes of entangled 
states: the GHZ and W Werner-type states. The former 
reads: 



Pg(0) = rpcHZ 



1-r 



-It 



8. 



(21) 



where p G nz=\ (|000> + |111)) «000| + (111|), r is the 
purity ranging from to 1, and 1% is the identity matrix of 
dimension 8. The explicit evaluation of the time-evolved 
state for both local and non-local system-environment 
coupling is reported in Appendix [X] The different fea- 
tures of the classical environmental noise of three qubits 
plays a key role into the dynamics of GHZ-type state. 
When the three qubits are coupled to different envi- 
ronments, the time-evolved density matrices takes an X 
shape as shown in Eq. (|A1[) . Instead when the interac- 
tion between the tripartite system and a common source 
of RTN is considered, all the elements of the three-qubit 
density matrix are non-vanishing after the initial time 
(see Eq. (|A3|) V Such a behavior differs from what was 
found in the two-qubit model [3l| where the time-evolved 
density matrices of any initial Werner state has an X 
shape regardless of local and non-local qubit-environment 
interaction [36j . 

In the other initial configuration, the three qubits are 
taken in the W Werner-type state 



1-r 

pw(0) = rp w H — 1 8 



(22) 



A. GHZ-type states 

1. Different Environments 

Here the effects of the dephasing on the state of 
Eq. (|2"Tj) are examined. When the qubits are coupled to 
different sources of RTN, the time-evolved state takes an 
X shape which allows one to calculate easily the tripar- 
tite negative and the expectation value of the EW given 
in Eq. (fT3"|). Their analytical forms are: 

max [0, ArG\(t) — (1 — r)] 



TT[p G (t)W W 2] = — 



ZrG\{t)- 



3-r 



(23) 



respectively. Such expressions indicate that both the 
tripartite negativity and the expectation value of the 
yVw2(t) can be expressed in terms of Specifically, 
we note that, at initial time for r ranging from 1/5 to 3/7 
the negativity is different from zero while the expectation 
value of Ww2 is zero or positive. This means that the 
EW does not detect the tripartite entanglement which is 
known to be present via the tripartite negativity mea- 
sure. This is in agreement with previous results in liter- 
ature 13]. In the limit of long times, G\{t) goes to zero, 
and, as a consequence Af^ vanishes and — (Ww2(t)} be- 
comes positive. In other words, the initial quantum cor- 
relations in the form of entanglement disappear in the 
limit of long times. 



withp„,=§ (|001) + |010) + 1 100)) ((001| + (010| + (100|). 
Unlike the GHZ Werner-type state, the time-evolved 
density matrix of the system, reported in Eqs. (|A5|) 
and (|A6[) takes the same form for local and non-local 
system-environment interaction. As it will be shown in 
the following, the peculiar dynamics exhibited by GHZ 
and W Werner-type states is crucial in suppressing or 
preserving the quantum correlations among the three 
qubits. 

It is worth noting that given the initial states of 
Eqs. pip and (|22|) and the form of the Hamiltonian writ- 
ten in Eq. (I15p , the time-evolved states are still invariant 
under the permutations of the three parties, regardless 
the local or non local nature of the system-environment 
coupling. This allows us to evaluate more simply both 
tripartite negativity and quantum discord. 



IV. RESULTS 

In this section, we evaluate the time evolution of the 
tripartite entanglement and quantum discord in the phys- 
ical model introduced in the previous section. The dy- 
namics of quantum correlations is investigated for the 
GHZ- and W-type states in the cases of local and non- 
local system-environment interaction. 



y/v=10 




FIG. 1. (Color online) Top panels: Negativity (left), oppo- 
site of the expectation value of the EW Ww2 (center), and 
tripartite quantum discord (right) as a function of jt and the 
purity r in the Markov regime with y/u= 10, when the three 
qubits, initially prepared in the state pa(0), are coupled to 
different sources of classical noise. Bottom panels: same as in 
the top panels in the non-Markov regime with 7/1/ = 0.1. 
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In Fig. [[J we report A/" (3 \ -(% 2 ), and X> (3) as a func- 
tion of the time and the initial purity of the state for two 
different values of j/v, namely 10 and 0.1, correspond- 
ing to the Markovian and non-Markovian dynamics of 
quantum correlations, respectively. In the first case, the 
RTN leads to a monotonic decay of entanglement and 
discord, while in the non-Markovian regime all quanti- 
ties are damped oscillating functions of time and display 
sudden death and revival phenomena. Now let us com- 
pare qualitatively the behavior of the three estimators 
of quantum correlations. Even if the tripartite negativ- 
ity and quantum discord present substantially the same 
qualitative behavior, nevertheless we find that V^) can 
take non- vanishing values in regions whereas A/"*- 3 -* is zero. 
As it can clearly be seen in the non-Markovian regime, 
the peaks of the tripartite negative found at a given value 
of jt decrease for smaller r's at higher rate than the cor- 
responding peaks of the tripartite discord, and, as con- 
sequence vanish at smaller values of the purity of initial 
state. Such a behavior results to be consistent with was 
found in bipartite systems of initially entangled qubits 
evolving under local decohering channels where separa- 
ble states can exhibit nonzero discord [l8|, [ljl |38[ . On the 
other hand, the entanglement and quantum correlations 
quantified in terms of Af^ and result to be higher, 
both in Markovian and non-Markovian regime, than the 
one detected by the tripartite entanglement witnesses. 



2. Common Environment 

Now let us evaluate the dynamics of the quantum cor- 
relations initially present in the state /?g(0) when the 
three qubits are coupled to a common source of RTN. 
The time-evolved tripartite negativity and expectation 
value of the witness Ww2 can be expressed as: 

A" (3) (i) = J max 0,2^2 (Gj(t) + 1) - (1 - r) 



Tr[ PG (t)W W 2] =-- 



\rGS)-(\-r 



(24) 



It is worth noting that, unlike the case of local system- 
environment interaction, here quantum correlations are 
not totally destroyed at long times. Indeed, the asymp- 
totic forms of Eqs. ([24| . holding for any degree of qubit- 
environment coupling, are given by 



lim JV {3) (t) = -ma,x [o, (2\/2 + l)r - 1 
- HmTr[p G (t)Ww2}= ( ^ - § 



(25) 



This means that the saturation values of A^ 3 -* and 
(/>gH , W2) do not depend upon the parameter 7 that is, 
in terms, upon the Markovianity of the regimes. For an 
initial state with purity r greater than 3/4, the entangle- 
ment quantified in terms of the negativity and that part 
of entanglement detectable by means of the witness Ww2 



do not vanish in the long time limit. The partial preser- 
vation of the entanglement can be ascribed to the indirect 
interaction among the qubits stemming from the coupling 
of the global system to a common noisy environment. 
Unlike the local system-environment interaction, here the 
environment is no longer only the source of decohering 
effects, but also represents a sort of interaction mediator 
between the subsystem. Such an interaction somehow 
hinders the destruction of quantum correlations. The 
survival of quantum correlations has already been ob- 
served for the bipartite entanglement of two-qubit sys- 
tems interacting with quantum environments 0, HI, l39| . 
When the purity ranges from 3/4 to 1/(2^2+ 1), the 
residual amount of tripartite entanglement at large times 
as quantified by A/"*- 3 ' is not detectable by means of the 
EWW K . 




FIG. 2. (Color online) Top panels: Negativity (left), opposite 
of the expectation value of EW Ww2 (center), and tripartite 
quantum discord (right) as a function of *yt and the purity r 
in the Markov regime with 7/^=10, when the three qubits, 
initially prepared in the state pghz, are coupled to a common 
source of classical noise. Bottom panels: same as in the top 
panels in the non-Markov regime with 7/1/ = 0.1. 



As shown in Fig. f2J also the tripartite quantum dis- 
cord can survive the decohering effects due to the RTN. 
In the Markovian regime, we find that all the estima- 
tors decay monotonically with time until they reach the 
corresponding saturation values depending upon r. In 
the non-Markovian dynamics, both A/"' 3 - 1 and exhibit 
damped oscillations but no sudden death of the correla- 
tions occurs for not-separable initial states. Indeed, the 
indirect interaction among the qubits not only counter- 
acts the total suppression of of negativity and discord at 
long times but also prevents the disappearance of quan- 
tum correlations at finite times. The survival of A/^ 3 -* 
and I?*- 3 ' in the long time limit represents the major dis- 
crepancy with what was found in the two-qubit form of 
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the model of Eqs. ([15]) and (JTSJ) 31]. There both bipar- 
tite entanglement and discord disappear at long times, 
regardless the local or non-local character of the interac- 
tion among the qubits and the environment. 

Finally, the time behavior exhibited by the detectable 
entanglement in the non-Markovian regime appears to be 
peculiar. Indeed, the expectation value of the EW Ww2 
takes positive or null values at finite times in correspon- 
dence of the non- vanishing minima of the tripartite neg- 
ativity. This implies that the suppression of the sudden 
death of entanglement cannot be detected while, as indi- 
cated by Eqs. ([25]) . the long-time entanglement protection 
can successfully be revealed. 



B. W-type states 

1. Different Environments 

The dynamics of the quantum correlations, initially 
present in the state pw(0) of Eq. (f22|) . is here analyzed 
for the case of qubit locally interacting with their envi- 
ronments. Unlike GHZ-type states discussed in the pre- 
vious section, the tripartite negativity can not be put in 
a compact analytical form, while the expectation value 
of the EW Wwi can be expressed as: 



y/v=10 




FIG. 3. (Color online) Top panels: Negativity (left), opposite 
of the expectation value of EW Wwi (center), and tripartite 
quantum discord (right) as a function of ft and the purity r 
in the Markov regime with j/v=10, when the three qubits, 
initially prepared in the state pw{0), are coupled to different 
sources of classical noise. Bottom panels: same as in the top 
panels in the non-Markovian regime with j/v = 0.1. 



Tr[p w (t)W 



wi\ 



-[r (relit) 



5Gl(t) +5G 2 (t) +4) - 13] . (26) 



At the initial time for r=l, the expectation value of Wwi 
is —1/3, and the tripartite negativity is equal to 0.94. 
The EW detects entanglement when r > 13/21, while 
W (3 ) is nonzero for r > 0.2096. 

Fig. [3] displays the tripartite negativity, the detectable 
entanglement and the discord as a function of time and 
initial purity of the state for two different values of 7/f, 
corresponding to Markov and non-Markov regime. In the 
first case, we find again the monotonic decay of A/"*- 3 - 1 and 
both vanishing at very long times. On the other 
hand the witness Wwi turns out to be unable to detect 
entanglement after a finite time. In the non-Markovian 
dynamics both the tripartite negativity and quantum dis- 
cord exhibit revivals and sudden death phenomena, while 
the entanglement detectable by means of Wwi still de- 
cays monotonically with time and vanishes at short times. 
Also for this qubits configuration, the state can be sepa- 
rable even with a non-zero discord. 



2. Common Environment 



EW Wwi takes the form 
Tv{pw{t)W W i] = 



32 



{9G 6 (t)+6G 4 (t) + 7G 2 (t)+6) 



13 
24 



(27) 



Unlike what was found in the previous section, at long 
times entanglement cannot be detected independently 
from the value of r. As it can be seen from Fig. 21 in 
the Markovian dynamics, the usual decay of quantum 
correlations is observed with tripartite negativity and 
quantum discord vanishing at long times. When the non 
Markovian dynamics is examined, A/"*- 3 - 1 and show 
damped oscillations without sudden death phenomena, 
even if they vanish in the long time limit (not displayed 
in the bottom panels of Fig. |4]). On the other hand, a re- 
vival and sudden death of the detectable entanglement are 
observed. These results together with the ones shown in 
the previous subsection clearly indicate that the quantum 
correlations initially present in the W-type states are less 
robust than those of the GHZ-type states. Indeed, both 
tripartite entanglement and quantum discord of the state 
pw(0) are completely suppressed not only in the case of 
subsystems each locally coupled to its environment but 
also for non-local system-environment coupling. 



Finally, we estimate the time evolution of the quantum 
correlations when the three qubits, initially set in the 
state pw{0), are subject to the decohering effects of a 
common source of RTN. The expectation value of the 



V. CONCLUSIONS 

Recently, a number of theoretical works examined the 
dynamics of non-classical correlations in tripartite sys- 
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FIG. 4. (Color online) Top panels: Negativity (left), opposite 
of the expectation value of EW Wwi (center), and tripartite 
quantum discord (right) as a function of yt and the purity r 
in the Markov regime with y/v= 10, when the three qubits, 
initially prepared in the state pw, are coupled to a common 
source of classical noise. Bottom panels: same as in the top 
panels in the non-Markovian regime with j/v = 0.1. 



tems coupled to external environments characterized by 
different memory properties [f3, [§-[3, While dif- 

ferent approaches have been adopted to estimate entan- 
glement [g, [{| [13, [H, [3 and Bell non- locality 11 1 , few 
studies focus on the time evolution of tripartite quantum 
correlations in terms of quantum discord |25j . 

In this paper, we address the dynamics of quantum cor- 
relations in a model consisting of three initially entangled 
qubits, not interacting among each other and coupled to 
different sources or to a common source of classical RTN. 
Specifically, two initial configurations have been consid- 
ered, namely GHZ and W Werner-type states. Unlike 
previous analyses [§-[l2j], here we have provided an es- 
timate of the tripartite quantum discord by using the 
approach recently introduced in Ref. [26j which allows 
one to quantify the amount of genuine tripartite quan- 
tum correlations of the system. These have then been 
compared with the entanglement evaluated by means of 
the tripartite negativity and with the detection ability of 
suitable EWs. 

Our results show that both entanglement and quan- 
tum discord are strongly affected not only by the initial 
configuration of the qubits but also by the local or non- 
local nature of the system-environment interaction. In 
particular, for a GHZ-type initial state the indirect in- 
teraction between the qubits due to their coupling to a 
common source of RTN allows for a long-time tripartite 
entanglement and quantum discord preservation, regard- 
less the Markov or non-Markov character of the envi- 
ronment itself. Furthermore, we find that the survived 



entanglement can be detected by means of the witnesses. 
On the other hand, quantum correlations are completely 
destroyed when the subsystems are coupled to different 
environments. The survival of entanglement and quan- 
tum discord and the disappearance of sudden death phe- 
nomena, which have already been observed in a number 
of bipartite systems interacting with quantum environ- 
ments [J, |g, |39| are here closely related to the tripartite 
nature of the system. Indeed, in the two-qubit model 
analogous to the one here investigated the number of 
quantum correlations vanish in the long time limit under 
any condition 13 II] . By using approaches already used 
elsewhere [29], l37l. l4dl. l4l| , it could certainly be of inter- 
est to extend our investigation to the case of 1// Q noises 
stemming from a collection of random telegraph sources 
with different switching rates. So, we could verify the 
survival of quantum correlations when a large number of 
decoherence channels is considered. 

In the other initial configuration, that is the qubits 
prepared in a W Werner-type state, no preservation of 
quantum correlations is found. Specifically, some of the 
standard features of Markov and non-Markov regime are 
observed. Indeed, the former exhibits a monotonic decay 
of quantum correlation while in the latter sudden death 
and revival phenomena occur only when the qubits are 
locally coupled to different sources of RTN. Anyway, both 
the entanglement and quantum discord of W Werner- 
type state turn out to be less robust than those of the 
GHZ-type state. 

Finally, the tripartite quantum discord deserves a brief 
comment. It almost shows the same qualitative behav- 
ior of the negativity for all the physical conditions ex- 
amined. Nevertheless, in agreement with what was pre- 
viously found in bipartite systems 0, [ll| HI, EH, the 
amount of tripartite quantum correlations quantified by 
discord can be different from zero for states without en- 
tanglement. From this point of view, our analysis seems 
to further validate the measure introduced in Ref. [26[ 
as a good estimator of the genuine tripartite quantum 
correlations which are not entanglement. 
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Appendix A: Evaluation of the time-evolved states 

Here, we give the explicit forms of the time-evolved 
state of the two different initial configurations of the 
qubits, namely GHZ and W Werner-type state, for the 
case of local and non-local system-environment interac- 
tion. 

In order to evaluate the dynamics of the system, first 
we calculate the evolution of the initial state for a given 
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choice of the noise parameter, as indicated in Eq. (fTTI) . 
Then, the obtained density matrix is averaged over noise 
(see Eq.([T9"])). For the input GHZ-type state p G (0) of 
Eq. (HH, we find that, in the case of local system- 
environment coupling, the time-evolved density matrix 
of the system takes the form 



with 



Pa(t) 



a(t) 




















a(t) 





9{t) 














m 











6{t) 








m 

















8{t) 


m 




















m 


9{t) 

















m 








9{t) 











m 














9(f) 





a(t) 




















a(t) 



a(t) = 
9(t) = 



= ; i 



3rG 2 2 (t)) 
l-rGa(t)), 



/8(t) = -(l-G^(t)), 
and a(t) = ^(l + 3Gi( t )). 

(A2) 



(Al) 



On the other hand, for the case of non local qubit- 
environment interaction the dynamics of the system ini- 
tially in /Og(0) results into: 



P N G L {t) = 



i + 3 M (i) 


-A(t) 


-AW 


-AW 


-AW 


-AW 


-AW 


§+3mW 


-A(t) 




AW 


A(t) 


AW 


AW 


AW 


-AW 


-A(t) 


A(*) 


I - mW 


A(i) 


AW 


AW 


AW 


-AW 


-A(t) 


A(*) 


A(t) 


I - mW 


AW 


AW 


AW 


-AW 


-A(t) 


A(t) 


A(t) 


A(*) 


I - mW 


AW 


AW 


-AW 


-A(t) 


A(t) 


A(t) 


AW 


AW 


1 - mW 


AW 


-AW 


-A(t) 


A(t) 


A(<) 


AW 


AW 


AW 


\ - mW 


-AW 


| + 3/i(*) 


-A(t) 


-AW 


-AW 


-AW 


-AW 


-AW 


| + 3 M W 



where 

^W = ^ (1 + G 4 W) , and AW = ^(1-G 4 W)- 

(A4) 

These results clearly show that the form of the time- 
evolved density matrix of three qubits initially in a GHZ- 
type state depend upon the local or non-local nature of 



the system-environment coupling. 

Now let us focus on the time evolution of the system 
initially prepared in the W-type state pwifi) given in 
Eq. ([^|) . When each qubit interacts locally with its en- 
vironment, the density matrix of the system at time t 
takes the form: 
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«(t) 



Pw(t) 



o 



«(*) 

12 



«(*) 
12 

<t) 
12 



r(t) | g £(«) 



24 



12 



12 ° 






x(t) 

12 



1*1 
12 



24 



12 






m 

12 



«(«) 
12 





24 T S 



y(«) 

12 

£(*) 
12 



£(£) 
12 

£(£) 
12 



24 






m 

12 



w(t) 
12 




12 



24 



y(*) 

12 



12 




fit) 

12 



<£(*) 

12 



£(*) 

24 



m 

12 
12 



m 

12 






x(t) 



whore 



J 



K(t)=r(l + G 2 (t)) (l-G 2 (t)), 

o 

r(t) = r(l + G 2 (t)) [(1 + G 2 (t)) 2 + 2(1 - G 2 (t)) 
m = r{l-G 2 {t))[2{l + G 2 (t)) 



1 



1(1 -r), X W = r(l-G 2 (t)) 2 (l + G 2 (t)), 

2 ], t (t)=r(l + G 2 (t))(l + G*(t)), 
p(t)=r(l-G 2 (t))(l + £?*(*)) 



When all the three qubits are coupled to a common 



i — 



source of RTN, the time-evolved state of the system can 
written as: 



5+A(t) 








H(t) 





S(i) 


S(t) 








I + T(t) 


n(*) 





n(t) 








r(i) 





fi(t) 


I + T(t) 





n(t) 








r(t) 


S(t) 













A(t) 


A(t) 








n(t) 


n(t) 





| + T(t) 








r(t) 


S(t) 








A(t) 





5+*(*) 


A(i) 





S(t) 








A(t) 





A(t) 


! + *(*) 








r(t) 


r(t) 





r(t) 








! + *(*) 



(A6) 



with 



A(t) = r 
T(i)=r 
r(t) = r 
A(i) = r 



1 3G 2 (t) 

16 32 

-1 7G 2 (t) 

48 96 



3G 4 (i) 3G 6 (t) 



1 

16 

5 

48 



3G 2 (t) 

32 
7G 2 (t) 

96 



16 

, G4(t) 
' 16 

G 4 (t) 

16 
G 4 ft) 

16 



+ 



32 
3G 6 (t) 

32 
3G 6 (f) 

32 
3G 6 (t) 
32 



S(t)=r 
0(f) = r 

*(*) = 
*(*) = 



1 

16 
5 

48 
-1 
48 " 

16 " 



3G 2 (t) G 4 (i) 3G 6 (t) 



32 
7G 2 (t) 

96 
7G 2 (t) 

96 
3G 2 ft) 

32 



+ 



16 
G 4 (t) 

16 
G 4 (t) 

16 
3G 4 (t) 

16 



32 
3G 6 (t) 

32 
3G 6 (t) 

32 

3G 6 (t) 
32 
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Unlike the GHZ Werner-type state, we note that in 
this case the time-evolved density matrix of the sys- 
tem takes the same form for local and non-local system- 
environment interaction. Indeed, in both the matrices of 



Eqs. (|A5[) and (| A6|) the diagonal 4x4 subblocks have an 
X shape, while in the antidiagonal ones the diagonal and 
antidiagonal elements are zero. 
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